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In this paper we report the first LDA+DMFT results (method combining Local Density Approx- 
imation with Dynamical Mean-Field Theory) for spectral properties of superconductor LaFePO. 
Calculated k-resolved spectral functions reproduce recent angle-resolved photoemission spectroscopy 
(ARPES) data [D. H. Lu et al, Nature 455, 81 (2008)]. Obtained effective electron mass enhance- 
ment values m* jrri ~ 1.9 - 2.2 are in good agreement with infrared and optical studies [M. M. Qazil- 
bash et al, Nature Phys. 5, 647 (2009)], de Haas-van Alphen, electrical resistivity, and electronic 
specific heat measurements results, that unambiguously evidence for moderate correlations strength 
in LaFePO. Similar values of m* /m were found in the other Fe-based superconductors with sub- 
stantially different superconducting transition temperatures. Thus, the dynamical correlation effects 
are essential in the Fe-based superconductors, but the strength of electronic correlations does not 
determine the value of superconducting transition temperature. 

PACS numbers: 71.27.+a, 71.10.-w, 79.60.-i 



I. INTRODUCTION. 

The discovery of superconductivity with the transi- 
tion temperature T c m 4 K in LaFePOi and T c s» 
26-55 K in ROi_ a F x FeAs (R = La, Sm) 2 has gener- 
ated great interest to the new class of Fe-based super- 
conductors.- While the microscopic mechanism of su- 
perconductivity in LaFePO is not yet clear^i 5 - its elec- 
tronic properties have been studied extensively. Vari- 
ous experiments revealed the presence of electronic cor- 
relations in LaFePO. Analyzing angle-resolved photoe- 
mission spectroscopy (ARPES) data, D. H. Lu et al. 
(Ref. 0) demonstrated that the DFT band structure 
should be renormalized by a factor of 2.2 to fit the exper- 
imental angle-resolved photoemission spectra. From in- 
frared and optical conductivity data the authors of Ref. ,7 
made a conclusion that the effective electron mass renor- 
malization is about 2 in LaFePO. Similarly, the elec- 
tron mass renormalization obtained from de Haas-van 
Alphen study is m*/m w 1.7 - 2.1 (Ref. 8) Compar- 
ison of the experimental electronic specific heat coeffi- 
cient 7„ = 10.1 mJ/mol K 2 for LaFePO 9 with the DFT- 
value 5.9 mJ/mol K 2 (Ref. llfj) gives the value of 1.7 for 
electron mass enhancement. Also the electrical resistiv- 
ity at low temperatures has T 2 -dependence^^ showing 
importance of correlation effects. However, the authors 
of Ref. [ill from analysis of x-ray absorption (XAS) and 
resonant inelastic x-ray scattering (RIXS) data for sev- 
eral iron pnictide compounds (SmOo.ssFeAs, BaFe2As2, 
LaFe2P2) arrived at a conclusion that correlations are 
not very strong here. 

So far only a few electronic structure calculations for 
LaFePO by the DFT-based first-principles methods with- 
out any account for electronic correlations have been 
reporte d 6 ' 10 i 12 ~ — thus making the problem of correlation 
effects study for this material very timely. 

The combination of Density Functional Theory (DFT) 



and Dynamical Mean-Field Theory (DMFT) called 
LDA+DMFT method 15 is presently recognized state-of- 
the-art many-particle method to study correlation effects 
in real compounds. LDA+DMFT calculations for the 
FeAs-based superconductors^— lead to diverse conclu- 
sions on the strength of electron correlations in these 
materials. The authors of Ref. [l^ have proposed an ex- 
tended classification scheme of the electronic correlation 
strength in the Fe-pnictides based on analysis of several 
relevant quantities: a ratio of the Coulomb parameter U 
and the band width W (U/W), quasiparticle mass en- 
hancement to* /to, k-resolved and k- integrated spectral 
functions A(Ts.,lj) and A(u>). Applying this scheme to 
LDA+DMFT results for BaFe2As2 they came to conclu- 
sion that this material should be regarded as a moderately 
correlated metal. In this work we report the results of 
LDA+DMFT study for electronic correlation effects in 
LaFePO. For this purpose we have calculated spectral 
functions A(k, w), effective electron mass enhancement 
to* /to and compare our results with the available mea- 
surements for LaFePO finding a very good agreement 
between calculated and experimental data. A moder- 
ate spectral functions renormalization corresponding to 
to* /to ks 2 was found in LaFePO similar to values ob- 
tained for BaFe 2 As2>i2' while superconducting transition 
temperature in those materials could be an order of mag- 
nitude different. 



II. METHOD. 

The LDA+DMFT scheme is constructed in the follow- 
ing way: First, a Hamiltonian Hlda is produced using 
converged LDA results for the system under investiga- 
tion, then the many-body Hamiltonian is set up, and fi- 
nally the corresponding self-consistent DMFT equations 
are solved. By projecting onto Wannier functions^ we 
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obtain an effective 22-band Hamiltonian which incorpo- 
rates five Fe d, three O p, and three P p orbitals per 
formula unit. In the present study we construct Wan- 
nier states for an energy window including both p and 
d bands. Thereby hybridization effects between p and 
d electrons were explicitly taken into account and eigen- 
values of the Wannier functions Hamiltonian H^da ex- 
actly correspond to the 22 Fe, O, and P bands from 
LDA. The LDA calculations were performed with the 
experimentally determined crystal structure^ using the 
Elk full-potential linearized augmented plane-wave (FP- 
LAPW) codei22 Parameters controlling the LAPW basis 
were kept to their default values. The calculated LDA 
band structure 6lda(^-) was found to be in good agree- 
ment with that of Lebegue et al. (Ref. fToh . 

The many-body Hamiltonian to be solved by DMFT 
has the form 

H = H LDA -H dc + - ^2 nf aa nfp a , , (1) 

where U£g' is the Coulomb interaction matrix, hf aa is 
the occupation number operator for the d electron with 
orbital a or /3 and spin indices a or a' in the i-th site. 
The term Hd c stands for the d-d interaction already ac- 
counted in LDA, so called double-counting correction. 
The double-counting has the form Hd c = U{ndm& — n)I 
where ridmft is the total self-consistent number of d elec- 
trons obtained within the LDA+DMFT and U is the av- 
erage Coulomb parameter for the d shell. 

The DMFT self-consistency equations were solved it- 
eratively for imaginary Matsubara frequencies. The aux- 
iliary impurity problem was solved by the hybridization 
function expansion Continuous-Time Quantum Monte- 
Carlo (CTQMC) method^ In the present implementa- 
tion of the CTQMC impurity solver the Coulomb inter- 
action is taken into account in density-density form. The 
elements of U°% matrix were parameterized by U and J 
according to procedure described in Ref. [24|. We used in- 
teraction parameters U — 3.1 eV and J = 1 eV similar to 
the values calculated by the constrained LDA method for 
Wannier functions^ in Fe-pnictides^ Calculations were 
performed in the paramagnetic state at the inverse tem- 
perature P = 1/T = 20 eV -1 . The real- axis self-energy 
needed to calculate spectral functions was obtained by 
the Pade approximant2& (see Appendix). 

III. RESULTS AND DISCUSSION. 

The orbitally resolved Fe 3d, O p and P p spectral func- 
tions computed within LDA and LDA+DMFT, respec- 
tively, are compared in Fig. [TJ Within the LDA all five 
Fe d orbitals form a common band in the energy range 
(—2.5, +2.0) eV relative to the Fermi level (band width 
W ~ 4.5 eV). There is a significant hybridization of the 
Fe 3d orbitals with the P p and O p orbitals, leading to 
appearance of Fe d states contribution in the energy inter- 
val (—5.5, —2.5) eV where the P p band is located. The 
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FIG. 1: (Color online) Orbitally resolved Fe 3d, O p and 
P p normalized spectral functions of LaFePO obtained within 
LDA+DMFT (upper panel) are compared with the LDA re- 
sults (lower panel). 



corresponding features of LDA+DMFT spectral func- 
tions (upper panel in Fig. [1]) in the energy area (—5.5, 
—2.5) eV should not be mistaken for Hubbard bands be- 
cause the same peaks are present in non-correlated LDA 
bands (lower panel in Fig. [IJ. Correlation effects do not 
result in Hubbard bands appearance but lead to signifi- 
cant renormalization of the spectral function around the 
Fermi energy: "compressing" of energy scale so that sep- 
aration between peaks of LDA+DMFT curves becomes 
s=s2 times smaller than in corresponding non-correlated 
spectra. 

It is instructive to plot energy dependence of real part 
of self-energy ReE(w) (see Fig. [5]). Peaks in spectral func- 
tion A(k, ui) are determined by the poles of (w — e(k) — 
S(w)) -1 function or the energy values uj = e(k) + Re£(cL>) 
(here e(k) is non-correlated band dispersion). In Fig. [5] 
together with ReS(w) a function uj + (Hd c )u is plotted 
as a stripe having the width of non-correlated band e(k). 
The peaks of spectral function A(\t, lo) correspond to en- 
ergy area where this stripe crosses ReE(w) curve. As one 
can see such crossing happens only once in the energy 
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interval around the Fermi level so that spectral functions 
will have only poles corresponding to quasiparticle bands 
and no Hubbard band poles will be observed. 
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FIG. 2: (Color online) The real part of the self-energy E(w) 
(black line) depicted together with uj + (H dc )a (light stripe, 
see text). 

A quantitative measure of the electron correlation 
strength is provided by the quasiparticle renormalization 
factor Z = (1 — f^U=o) _1 which gives an effective mass 
enhancement m* jvn = . In general, the self-energy is 
a matrix, leading to different effective masses for different 
bands. The calculated m* jm values for every (i-orbital 
are presented in Table 1. The d x 2_ y 2 orbital has the 
smallest effective mass renormalization m*/m = 1.942. 
The other d orbitals have approximately the same value 
m*/m « 2.2. 

TABLE I: Effective mass renormalization m*/m of quasipar- 
ticles in LaFePO for different orbitals of the Fe d shell from 
the LDA+DMFT calculation. 



Orbitals 


dxy dy ZtX z d^ z 2 _ r 2 d x 2_y2 


m*/m 


2.189 2.152 2.193 1.942 



The calculated effective mass enhancement 
m*/m « 1.9 - 2.2 in LaFePO agrees very well 
with the de Haas-van Alphen experiments^ where it was 
found to range from 1.7 to 2.1, and with the estimations 
of the effective mass renormalization of a factor of 
2 from optical conductivity data 7 and specific heat 



measurements. 

We now calculate the k-resolved spectral function 

A(k,u) = -Im-Tr[(u> + n)I -h^-tiw)}" 1 . (2) 

Here = Hlda — Hd c is the 22x22 Hamiltonian matrix 
on a mesh of k-points and \x is the self-consistently deter- 
mined chemical potential. In Fig. [3] we compare our re- 
sults with ARPES data of Lu et al. (Ref. M) . Both theory 
and experiment show dispersive bands crossing the Fermi 
level near the T and M points. In addition, two bands 
can be seen at -0.2 eV and in the region from -0.3 to -0.4 
eV near the T point. The calculated shape and size of the 
hole and electron pockets centered at the T and M points, 
respectively, are in good agreement with the ARPES, see 
Fig. [3] (lower panel), and de Haas-van Alphen^ data. 
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FIG. 3: (Color online) The k-resolved total spectral function 
A(k, lu) of LaFePO along the F - X - T and F - M lines 
in the Brillouin zone is depicted as a contour plot. Upper 
panel: The LDA+DMFT spectral function. Lower panel: 
The corresponding experimental ARPES intensity map of 
Lu et al. (Ref.H). 

The correlated band structure cdmftO*.) is also shown 
in Fig. [3] (upper panel). Near the Fermi energy, i.e., 
in the energy range from -0.2 eV to zero where quasi- 
particles are well defined (as expressed by a linear be- 
havior of ReE(w), see Fig. [2]), this dispersion is very 
well represented by the scaling relation e d mft (k) = 
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£iDA(k)/(m*/m), with m*/m taken as the computed 
mass enhancement from Table 1. 

The present results mean that the band structure in 
LaFePO is renormalized by the correlations and con- 
sists of the quasi-two-dimensional Fermi surface sheets 
observed in the experiments. At the same time, there is 
no substantial spectral weight transfer from quasiparti- 
cle bands near the Fermi energy to Hubbard bands and 
the system is far away from metal-insulator transition. 
Such spectral function behavior does not allow to classify 
LaFePO either as strongly or weakly correlated material. 
According to classification proposed in Ref. [l^ this ma- 
terial can be regarded as a moderately correlated metal, 
like BaFe 2 As 2 J£ 

In other iron pnictide materials the electron mass en- 
hancement was also reported to be close to value of 
«2. For example, from the de Haas-van Alphen ex- 
periments m* /m — 1.13 - 3.41 (Ref. 27) was found for 
SrFe2P2, isostructural analogue of the superconducting 
compounds Sri_ x K x Fe2As2 (T c = 37 K). For the series 
of isostructural compounds Ba(Fei_ x Co :c )2As2, x = 
0.3 the authors of Ref. [28| evaluated m*/m = 2-4 from 
the ARPES data. Also from the analysis of ARPES 
an electron mass renormalization by 2.7 was evaluated 
in Bao.6Ko.4Fe2As2 (T c = 37 K) ?9 as comparing with 
the DFT-band structure. Further, the DFT-calculated 
plasma frequencies are by a factor of 1.5 to 2 larger than 
the experimental values^ showing the electron mass en- 
hancement of the same strength for LaFePO, LaFeAsO, 
SrFe 2 As 2 , BaFe 2 As 2 , K .45Bao.55Fe 2 As2 (T c « 30 K), 
and LaOo.gFo.iFeAs (T c » 26 K). From the reported 
data, it follows that the electronic correlation strength in 
the FeAs-based superconductors with the superconduct- 
ing transition temperatures up to 37 K and their parent 
compounds is the same as in superconductor LaFePO 
with T c w 4 K. Thus, the strength of electronic corre- 
lations in the Fe-based superconductors seems to be not 
intrinsically connected with the superconducting transi- 
tion temperature. 



IV. CONCLUSION. 

By employing the LDA+DMFT method we have calcu- 
lated the spectral functions and single-particle k-resolved 
spectrum of superconductor LaFePO for the first time. 
Very good agreement with the ARPES data was found. 
In the spectral functions we observed no substantial spec- 
tral weight transfer. The obtained effective electron 
mass enhancement values m*/m w 1.9 - 2.2 are in good 
agreement with infrared and optical studies, de Haas- 
van Alphen, and specific heat results. The electronic 
correlation strength in LaFePO with small value of su- 
perconducting temperature 4 K is similar to the other 
Fe-pnictide superconductors with the transition temper- 
atures up to 37 K. 
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VI. APPENDIX: SELF-ENERGY ON THE REAL 

AXIS 

In the DMFT method temperature (Matsubara) Green 
functions formalism is used with arguments in the form 
of imaginary time r or corresponding Matsubara imagi- 
nary energies iui n = i(2n + 1)tt//L In order to calculate 
a spectral function one needs to have self-energy as a 
function of real energy that means to perform an- 

alytical continuation of the function Y>(iuj n ) to real axis. 
One of the usual algorithms for analytical continuation 
is Pade approximantfSS, and it was successfully used in 
earlier LDA+DMFT calculations where effective impu- 
rity problem was solved by Iterative Perturbation The- 
ory (IPT) methodi^I However, when impurity problem 
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FIG. 4: (Color online) The imaginary part of the Fe d yz self- 
energy £(iw n ) calculated within LDA+DMFT (black circles) 
is compared with the corresponding Pade approximant (red 
curve) . 
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FIG. 5: (Color online) Orbitally resolved Fe 3 d spectral func- 
tions of LaFePO from the Maximum Entropy method (green 
shaded areas) and the real-axis self-energy £(w) obtained with 
the use of Pade approximation (blue curves). 



is solved by stochastic Quantum Monte Carlo (QMC) 
method numerical noise appears in calculated E(iw„), 
see Fig. S) Attempts to apply the Pade approximant 
method to such noisy data result in completely wrong 
widely oscillating real axis function. 

In order to solve this problem Maximum Entropy 
method (MEM) method-^ was proposed. In this method 
spectral function A(u) corresponding to imaginary time 
Green function G(r) from QMC calculation is found as a 
best approximation to solution of the integral equation: 



G(r) 



duj- 



A(w) , 



(3) 



with the condition of maximization of effective entropy 
functional that gives a smooth spectral function. The re- 
sulting spectral function A(u) is identified then with the 
k-integrated analogue of Eq. [2] that gives equations for 
unknown self-energy E(u>). For many-orbital case that 
gives a set of equations with the corresponding number 
of unknown variables Ej(o;). Solution of such a set of 
equations can be a rather difficult problem. In addition 
to that the MEM method smears out all high-energy fea- 
tures in A{uj) due to the factor e - ^" in the kernel of 
integral equation ([3]). 

In the present work we have used a modified version 
of the Pade approximant method. To make the analyti- 
cal continuation procedure of the noisy self-energy E(iu>) 
numerically stable, we construct the approximant using 
only those frequencies values where the self-energy is a 
smooth function. Practically, that means to use a few 
first iuj n at the lowest frequencies and the data at the 
large frequencies where E(iw) approaches asymptotic be- 
havior. In the result we obtain a smooth function that 
has correct analytical behavior at two limits: small en- 
ergies uj — > and large energies where it obeys known 
asymptotic. 

In Fig.|4]we compare the Fe d yz self-energy E(zw n ) ob- 
tained within Pade approximation with the correspond- 
ing numerical data from QMC solution of DMFT equa- 
tions. The approximant has accurate derivatives in vicin- 
ity of zeroth Matsubara frequency, that guaranties cor- 
rect analytical properties near the Fermi level and correct 
asymptotic behavior. At the same time it approximates 
the noisy region with a smooth curve. 

In Fig. [5] the Fe 3d spectral functions obtained with 
Pade real-axis self-energy are compared with the MEM 
curves. The results for energies near the Fermi level are 
in very good agreement with each other. However, going 
to the higher and lower energies MEM curve very soon 
becomes smeared and nearly featureless while the curve 
obtained with Pade real-axis self-energy has much better 
resolved peaks and shoulders. For example, all the peaks 
in the energy region (—6, —2) eV are completely missed 
in MEM curve. This is due to exponential nature of the 
MEM kernel that suppresses all high-energy features. 



1 Y. Kamihara et al, J. Am. Chem. Soc. 128, 10012 (2006). 

2 Y. Kamihara et al, J. Am. Chem. Soc. 130, 3296 (2008); 
Z.-A. Ren et al, Chinese Phys. Lett. 25, 2215 (2008). 

3 Yu. A. Izyumov and E. Z. Kurmaev, Physics-Uspekhi 51, 
23 (2008). 

4 K. Kuroki et al, Phys. Rev. B 79, 224511 (2009); 
J. D. Fletcher et al, Phys. Rev. Lett. 102, 147001 (2009). 

5 M. Yamashita et al, Phys. Rev. B 80, 220509(R) (2009). 

6 D. H. Lu et al, Nature (London) 455, 81 (2008); 
D. H. Lu et al, Physica C 469, 452 (2009). 

7 M. M. Qazilbash et al, Nature Phys. 5, 647 (2009). 

8 A. I. Coldea et al, Phys. Rev. Lett. 101, 216402 (2008). 

9 S. Suzuki et al, J. Phys. Soc. Jpn. 78, 114712 (2009). 



S. Lebegue, Phys. Rev. B 75, 035110 (2007). 

1 W. L. Yang et al, Phys. Rev. B 80, 014508 (2009). 

2 V. Vildosola et ai, Phys. Rev. B 78, 064518 (2008). 

3 R. Che et al, Phys. Rev. B 77, 184518 (2008). 

4 Y. Kamihara et al, Phys. Rev. B 77, 214515 (2008). 

5 G. Kotliar et al, Rev. Mod. Phys. 78, 865 (2006); K. Held 
et al, Psi-k Newsletter 56, 65 (2003), reprinted in Phys. 
Status Solidi B 243, 2599 (2006). 

6 K. Haule, J. H. Shim, and G. Kotliar, Phys. Rev. Lett. 
100, 226402 (2008). 

7 M. Aichhorn et al, Phys. Rev. B 80, 085101 (2009). 

8 V. I. Anisimov et al, J. Phys.: Condens. Matter 21, 075602 
(2009); A. O. Shorikov et ai, JETP 108, 121 (2009); 



6 



V. I. Anisimov et al, Physica C 469, 442 (2009). 

19 S. L. Skornyakov et al, Phys. Rev. B 80, 092501 (2009). 

20 V. I. Anisimov et al, Phys. Rev. B 71, 125119 (2005). 

21 B. I. Zimmer et al, J. Alloy s Compd. 229, 238 (1995). 

22 http://elk.sourceforge.net/ 

23 P. Werner et al, Phys. Rev. Lett. 97, 076405 (2006). 

24 A. I. Liechtenstein, V. I. Anisimov, and J. Zaanen, Phys. 
Rev. B 52, R5467 (1995). 

25 Dm. Korotin et al, Euro. Phys. J. B 65, 1434 (2008). 

26 H. J. Vidberg and J. W. Serene, J. Low Temp. Phys. 29, 
179 (1977). 



27 J. G. Analytis et al, Phys. Rev. Lett. 103, 076401 (2009). 

28 V. Brouet et al, Phys. Rev. B 80, 165115 (2009). 

29 M. Yi et al, Phys. Rev. B 80, 024515 (2009). 

30 S. L. Drechsler et al, arXiv:cond-mat/0904.0827. 

31 V. I. Anisimov, A. I. Poteryaev, M. A. Korotin, 
A. O. Anokhin, G. Kotliar, J. Phys.: Condens. Matter 
9, 7359 (1997). 

32 M. Jarrell and J. E. Gubernatis, Physics Reports 269, 133 
(1996). 



